Scattering by an oscillating barrier: quantum, classical, and semiclassical comparison 
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We present a detailed study of scattering by an ampHtude-modulated potential barrier using 
three distinct physical frameworks; quantum, classical, and semiclassical. Classical physics gives 
bounds on the energy and momentum of the scattered particle, while also providing the foundation 
for semiclassical theory. We use the semiclassical approach to selectively add quantum-mechanical 
effects such as interference and diffraction. We find good agreement between the quantum and 
semiclassical momentum distributions. Our methods and results can be used to understand quantum 
and classical aspects of transport mechanisms involving time-varying potentials, such as quantum 
pumping. 
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I. INTRODUCTION 

Scattering dynamics involving periodic time-varying 
potentials is of fundamental importance to quantum 
transport physics and related applications in mesoscopic 
condensed matter physics. The quantum-mechanical 
treatment of an oscillating barrier was first studied by 
Biittiker and Landauer in order to understand electron 
tunneling times [1], and their work built on previous work 
on photon-assisted tunneling in superconducting diode 
junctions [2j. Since then several workers have developed 
theoretical tools for treating time- varying barrier or well 
potentials, for studying photon-assisted tunneling |3| , 
quantum pumping [5] , and electron scattering by intense 
laser-driven potentials j6]. These systems can display 
rich quantum and classical dynamics that include chaotic 
scattering and chaos-assisted tunneling iTrtllJ. dvnamical 
localization [l^ . and quantum interference |13| . 

Scattering by an amplitude-modulated potential bar- 
rier is of fundamental interest on its own, and it is also 
a building block for the more complex time-dependent 
potentials used in quantum pumping fT^ - [T6j . For ex- 
ample, the turnstile pump employs two potential bar- 
riers whose amplitudes oscillate tt/2 out of phase from 
each other. Despite its technological promise of generat- 
ing highly controlled and reversible currents at the single 
electron level [17|, quantum pumpin g in normal meso- 
scopic conductors remains elusive [isl. Il9l| (though it has 
been recently observed in a hybrid superconducting sys- 
tem [13). 

Experimental systems based on ultracold atoms offer 
the possibility of conducting precision tests of quantum 
pumping theories, while avoiding the capacitive coupling 
and rectification effects that have plagued attempted 
solid state implementations [19]. Furthermore, the use of 
ultracold atomic gases allows control over the momenta 
of the pumped particles and the coherence of the gas, 
and it permits precision imaging and velocity measure- 
ments, as well as the choice between Bose-Einstein and 



Fermi-Dirac statistics. 

In this paper, we study classical, semiclassical, and 
quantum dynamics of one-dimensional scattering by an 
amplitude-modulated Gaussian barrier. Motivated by 
possible experimental implementations with ultracold 
atoms, our main theoretical results are based on cal- 
culations of the scattered momentum distribution for 
atomic wavepackets of well-defined incident velocity, such 
as propagating Bose-Einstein condensates (BEG). By em- 
ploying a semiclassical formalism, we start with the clas- 
sical dynamics and selectively turn on quantum processes 
such as interference and diffraction. Our main results 
can be summarized as follows: i) Glassical physics gives 
bounds on the range of scattered momentum states; ii) 
Semiclassical and full quantum calculations predict simi- 
lar final momentum distributions; (iii) The heights of Flo- 
quet peaks, which are not easily predicted by quantum 
calculations, are explained quantitatively by the semi- 
classical method. Interestingly, the physical pictures for 
the scattering process are quite different for the semiclas- 
sical and quantum methods. The semiclassical approach 
interprets the discrete final momentum values as inter- 
cycle interference over multiple barrier oscillations, but 
with the relative amplitudes of these states determined 
by intra-cycle interference. In contrast, from the Floquet 
perspective of full quantum theory, the final momentum 
states can be viewed as sidebands of the initial momen- 
tum state. 

The paper is structured as follows: We present our 
model in Sec. |lTl and in Sec. IIIII display results of quan- 
tum and classical calculations for this model. Sec. IIVI 
explains the algorithm used for the semiclassical calcula- 
tion, and Sec. IVl compares and discusses the semiclassical 
and full quantum methods. In Sec. IVIl we show how the 
model and results of this paper can be tested experi- 
mentally with ultracold atoms. Sec. IVIl] summarizes our 
main results. Appendices |A] and |B] fill in the details of 
the semiclassical algorithm, and Appendix[C]explains the 
range of scattered momenta based on a simpler potential. 
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II. MODEL 

Our model is motivated by recent proposals [111, 
to simulate mesoscopic transport processes by studying 
ultracold atomic wavepackets propagating in quasi-one- 
dimensional waveguides that scatter from well-defined, 
localized potentials. A laser beam, blue-detuned from 
an atomic resonance, and tightly focused at the center 
of the wave guide, can create a potential barrier with a 
Gaussian profile, its width determined by the laser spot 
size and its amplitude by the intensity of the laser. 

We choose a ID Gaussian barrier, centered at the ori- 
gin, whose amplitude is modulated sinusoidally at fre- 
quency oj, with potential energy U{x,t) given by 



U{x, t) ^Uo{l + A sin(wt + (/)))e 



(1) 



Uq is the average amplitude of the barrier, A is the rel- 
ative modulation amplitude, a is the standard deviation 
width of the barrier, and (p is the phase of the modu- 
lation. The Hamiltonian describing particle motion and 
scattering from this potential is 



H 



P 

2m 



U(x,t). 



(2) 



We use wavepackets with initial momentum > 0, 
centered at a point x far to the left of the barrier, and 
whose position-space wavefunction is given by 



= 0) ^F(x)e 



iPfyX 



(3) 



where F{x) is the envelope of the wavepacket and is typ- 
ically a Gaussian of width /3, 



F{x)^Fg{x) = 



1 



(2^) 



1/4^ 



(4) 



Alternatively, the envelope may have a Thomas-Fermi 
distribution of radius /3, such that 



F{x) = Ftf {x) = 



ip'^~{x-xy ,\x-x\<p 

,\x-x\>f3 



(5) 



The Thomas-Fermi and Gaussian envelopes are typical 
of BEG wavefunctions in strongly interacting and non- 
interacting limits, respectively. Unless otherwise noted, 
we employ wavepackets that are much wider than the 
barrier width (/3 ^ cr), with packet width /? sufficiently 
large such that /3 ^ 2'KpQ/mLjj, ensuring that many bar- 
rier oscillations occur while the packet interacts with the 
barrier. 

In the rest of the paper, unless otherwise mentioned, 
we use Uq — m = h = 1^ A = 0.5, a = 10, and /3 = 300. 
The values of the incident momentum are in the range 
Pq ~ 1 — 2, the oscillation frequency uj ~ — 0.2, and 
in most cases the phase, (p, is set equal to 0. In the case 
of a Gaussian packet, we select x — —1500 to ensure 
separation of the initial packet from the barrier. 
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FIG. 1: (Color Online) Snapshots from a typical quantum- 
mechanical calculation showing a Thomas- Fermi (Eq. ©) 
wavepacket (left axis; solid red line) scattering off a Gaus- 
sian barrier (right axis; dotted green line). The amplitude 
of the barrier varies in time according to Eq. ^ with 
Uo = 1, A = 0.9. The barrier width (cr = 10) here is typi- 
cal in our simulations, but the packet width (/? = 40) is much 
less (to show more details) than used elsewhere (/3 = 300) in 
the paper. 



The choice of a theoretical unit convention based on 
h = 1 and m = 1 is equivalent to selecting an arbitrary 
time unit i„ and a related length unit lu ~ \/htu/m, with 
h — 1.054 X 10"'^'* J-s. The corresponding energy unit is 
Eu = f^/tm while the mass unit is that of the particle, 
TO,, = m. 



III. QUANTUM AND CLASSICAL 
CALCULATIONS 



1. Quantum Description 

We consider both quantum-mechanical and classical 
descriptions of the scattering process. This dual frame- 
work allows us to distinguish the classical and quantum 
nature of a variety of scattering features. 

Our quantum-mechanical approach is based on propa- 
gating the wavepacket with the Schrodinger equation: 



2to 



dl^ + U{x,t)m ^ -ihdt^ 



(6) 



via a split-step operator method |23| that incorporates 
the time- variation of the scattering potential U{x, t). The 
numerical calculation is done using a Fast-Fourier Trans- 
form (FFT) in a parallefized routine in FORTRAN. With 
periodic boundary conditions implicit in the FFT, the 
spatial range R (typically ~ 8000 in dimensionless units) 
is chosen sufficiently large to allow the entire wavepacket 
to interact with the barrier at R/2 without significant 
wraparound. The spatial grid density and the time step 
for propagation are both taken to be of the order of 0. 1 in 
dimensionless units. The resulting momentum grid den- 
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sity 2tt/R ~ 10~^ is more than sufficient to resolve the 
narrowest momentum space features that we encounter. 

Figure [1] shows a quantum calculation of a Thomas- 
Fermi wavepacket in position space at four separate times 
as it scatters from an amplitude-modulated Gaussian 
barrier. In order to show more details of the scattering, 
the packet width shown in this figure is intentionally more 
narrow than that used in the rest of the paper. The re- 
sulting transmitted and reflected wavepackets show con- 
siderable structure, but with no clear pattern, except for 
some residual spatial oscillation suggesting some type of 
interference effect. While examining the scattering pro- 
cess in position space does not yield any simple clues re- 
garding its dynamics, the momentum-space picture offers 
significantly more insight into the relevant physics. 

To obtain the wavefunction in momentum space, at a 
chosen large time, t = tf, after the packet has moved 
away from the potential barrier, we compute the Fourier 
transform of 'i'{x,tf) : 



1 



''^ {x,tf)dx. 



(7) 



We also compute the corresponding final-momentum 
probability density, 



P^{pf) = \^{Vf,tf) 



(8) 



Here, p/ is used to indicate momentum at the chosen 
final time. Also, we note that for sufficiently large times, 
such that the packet has moved far from the barrier, the 
final momentum distribution is constant in time, while 
the momentum-space wavefunction is not. 

A time-periodic potential produces energy and momen- 
tum sidebands to the incident carrier momentum state, 
which can be described by Floquet theory, the temporal 
analog of Bloch's theorem. In our model, a wavepacket 
is incident on the barrier with fixed group momentum 
po and associated kinetic energy = Pq/(2to). Since 
we use spatially-broad packets, the incident packet has a 
very narrow momentum spread. The interaction of the 
incident wavepacket with the amplitude-modulated bar- 
rier produces a series of discrete momentum states sep- 
arated in energy by hw. The allowed final-momentum 
states must obey the equation 



Pf{n) — ±y^2m (Eq + nhbj) 



(9) 



where n is any integer satisfying n > —Eo/hiu, and with 
and (— ) corresponding to transmission and reflec- 
tion, respectively. 

Fig. [2] shows the momentum-space distribution of 
the reflected and transmitted wavepackets after scat- 
tering from the amplitude-modulated barrier. The re- 
sults of the full quantum calculation show the regular 
"comb" of discrete momentum states consistent with 
Eq. Q. Fig. [2] also plots the classical momentum- 
space distribution for a Gaussian ensemble of particles 
with the same initial momentum spread as the initial 
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FIG. 2: (Color Online) Quantum (sharply peaked curves, blue 
online) and classical (green online) momentum distribution 
for fixed lo = 0.1, Uo ^ I, A = 0.5, cr = 10, /3 = 300, but 
different incident packet velocities, po- The classical distribu- 
tions were obtained via the histogram method, and statistics 
account for the fluctuations seen in the curves, (a) The re- 
flected and (b) transmitted parts for po = 1.4142; (c) reflected 
part for p = 1.0, when transmission is negligible; (d) trans- 
mitted part Po = 1-8, when reflection is negligible. 



quantum wavepacket (see next subsection for details). 
The classically-allowed bounds for the final momentum 
roughly constrain the Floquet comb on both refiection 
and transmission, though we find that the comb often ex- 
tends slightly past the classically-allowed bounds. How- 
ever, the amplitude of the teeth of the comb do not ap- 
pear to have any obvious pattern, and only loosely follow 
the strength of the classical final- momentum distribution. 

The semiclassical approach presented in Sec. IIVI and 
the Appendices will provide an alternative explanation 
for the positions of the teeth of the Floquet comb in 
terms of inter-cycle interference, and will provide an ex- 
planation for the relative amplitudes of the comb teeth 
in terms of infra-cycle interference. 



2. Classical Description 

The classical description of the scattering dynamics 
computes trajectories based on the Hamiltonian of Eq. 

In the static limit, particles of incident energy above 
Uq are transmitted, and those below are reflected. In 
contrast, scattering from an oscillating barrier leads to 
significant changes in the particle momentum distribu- 
tion, as particles gain or lose energy with the rise and 
fall of the potential. The final outcome depends on the 
phase of the oscillation as the particle encounters the 
barrier, and must generally be computed numerically. 

Our quantum and semiclassical calculations diminish 
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FIG. 3: (Color Online) (a) Final momentum vs. initial posi- 
tion for A = 0.5, Lj = 0.1, pq — 1.8. Capital letters correspond 
to different momentum regions (separated by solid vertical 
lines; see App. [Aj. At a selected pf, marked by the dashed 
line, paths arrive after beginning at many different xo; those 
points are labeled by Greek letters. Each lies on a branch 
of the multivalued function xo{pf,tf), and each branch is la- 
beled by a Roman letter, (b) Final-momentum distributions 
calculated quantum-mechanically and classically. The classi- 
cal calculations show a histogram (solid line, red online) and 
PciPf) from Eq. pip (dashed curve, black online). Fluctua- 
tions in the histogram arise for statistical reasons. 



the role of the phase of the barrier oscillation by study- 
ing Heisenberg-limited wavepackets with a large position 
spread and a well defined momentum, so that many bar- 
rier oscillations occur while the wavepacket is interacting 
with it. We mimic such wavepackets in our classical ap- 
proach by employing ensembles of particles with initial 
conditions whose position and momentum distributions, 
(x) and Pq (p), match those of the quantum distribu- 
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FIG. 4: (Color Online) Momentum distributions for fixed ve- 
locity of incident packet but for different values of cj. The 
correlation between classical and quantum distributions re- 
flected in Fig. [5] is seen. Comparison of quantum (blue and 
above axis) and classical (green and below axis) momentum 
distributions for po = 1.0 with oj = 0, 0.00263, 0.2, 0.8. Quan- 
tum and classical results are correlated for low and high values 
of uj, with significant differences appearing at intermediate 
values. These correlations are quantified in Fig. [5] 



tions: 



(10a) 
(10b) 



Generally, our initial momentum distributions are suf- 
ficiently narrow that classical particles can begin with a 
fixed initial momentum, distributed along a line segment 
that substantially covers the width of the initial wave 
packet, with statistical weights P^{x). 

The distribution PciPf) of final momenta pf can 
be obtained by numerically integrating trajectories and 
grouping them in bins of final momentum to plot a his- 
togram, as shown in Fig. [5] and Fig. [S^b). Alterna- 
tively, we can compute trajectories numerically to ob- 
tain the final momentum as a function of initial posi- 
tion xq and final time tf^ pf — p{x(),tf), as shown in 
Fig. ^a) . We note that due to the periodicity of the bar- 
rier amplitude, is a continuous periodic function of xq, 
with period 2T:pQ/uj'm. Any such periodic function has 
a maximum and minimum, which define the classically- 
allowed range oipf, as shown in Fig. [31 Furthermore, this 
periodicity means that many initial positions XQ{pf,tf) 
contribute to the final momentum distribution P^{pf). 
Each XQ{pf,tf) contributes to PciPf) ^ term propor- 
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FIG. 6: (Color Online) Phase sensitivity for very low omega 
values: even for the same uj = 0.0125 the momentum distri- 
bution changes with phase of the barrier oscilation: (j> = 
(blue, above axis) and = tt (green, below axis) with (a) 
showing transmitted fraction and (b) reflected fraction. The 
incoming wave packet was Gaussian-shaped, with po ~ 1.7 
and j3 — 300, and the barrier parameters were a = 10 and 
A = 0.5. 



FIG. 5: (Color Online) Correlation coefficient XQC of the nor- 
malized classical and quantum momentum densities (defined 
in Eq. (|12|l 'l plotted as a function of the barrier oscillation fre- 
quency LJ. The three different incident momentum represent: 
(i) Primarily transmitting (po = 1.8), (ii) primarily reflecting 
{po = 1.0) and (iii) a transition regime (po = 1.4142) of partial 
reflection and partial transmission. All use Gaussian packets, 
except that for po — 1.0. The result of using a Thomas-Fermi 
packet is also plotted, showing that the choice of packet shape 
is not crucial if the packets are sufficiently broad. 



tional to \dxl/dpf \ = \dp{xo,tf)/dxo\ 



'xo=3:lipf,tf)' 



iPf)-J2Pc{4iPf^tf))\d4/dpf 



(11) 



Figure [3] shows the final classical momentum distribu- 
tion PciPf) computed by both the histogram method 
(solid line, red online) and according to Eq. (jTl]) (dashed 
curve, black online), as well as the final quantum mo- 
mentum distribution. The maximum and minimum oipf 
define the classically-allowed region, with dpf/dxo going 
to zero at these locations, and its reciprocal in Eq. ((TT]) 
tending to infinity ^24.] . 

When we compare the quantum calculation to this 
classical calculation (Figs. [5] and [3Ib)) we see that the 
boundaries of the classically-allowed region accurately 
define the region of momentum space in which Floquet 
peaks are large. Small peaks also appear outside but 
close to the classically-allowed region. As we show in the 
semiclassical treatment of Sec. IIVI these are the result 
of momentum-space tunneling (or diffraction) into the 
classically-forbidden region. 

We also find that the barrier oscillation frequency oj, 
an easily variable experimental parameter, can be used 
to control the concurrence of the classical and quantum 
calculations, with good agreement in the limits of very 
high and low frequencies. For a static barrier, or for 
extremely low frequencies, momentum conservation in 
classical and quantum theories ensures agreement. As 



the frequency is increased, keeping the initial packet un- 
changed, the agreement gets poorer (Figs.|3]Jb) andjU^b)). 
The classical momentum distribution broadens, and the 
quantum distribution acquires a "comb" structure since 
Floquet peaks begin to resolve as their separations be- 
come greater than their widths (which depend inversely 
on the width of the initial packet in position space) . This 
is the range of particular interest in this paper. At very 
high frequencies, the incident particles cannot respond 
fast enough to the modulation of the barrier, and so they 
effectively interact with the time-average of the potential. 
The classically-allowed region narrows, while in the Flo- 
quet picture, the spacing between the Floquet peaks in- 
creases (Fig.m^c)). When there is only one non- negligible 
Floquet peak remaining, it coincides with the classically- 
allowed region, resulting again in good agreement be- 
tween the two methods (Fig. UJ^d)). 

In order to quantify the comparison of the final quan- 
tum momentum distribution, Eq. ([5]) , with the final clas- 
sical momentum distribution, Eq. (jlip . we define a kind 
of final momentum-density correlation coefficient. 



Xqc — 



JdpfP^iPf)PS{pf) 



jdpf P^ipf) jdpf PEiPf) 



(12) 



This correlation coefficient is plotted in Fig. [SJ which 
confirms the behavior indicated above, wherein the quan- 
tum and the classical distributions are in close agreement 
at low and at high frequencies, but not at intermediate 
frequencies. 

At very low frequencies, the quantum momentum dis- 
tribution depends upon the initial phase (j) of the po- 
tential (see Eq. ([1])), as we show in Fig. [HI At such low 
frequencies, the incident packet interacts with the barrier 
over only a fraction of a cycle, thus experiencing a barrier 
amplitude that is strongly dependent on the oscillation 
phase. 

To summarize, we see that classical calculations de- 
scribe the range of momenta over which Floquet peaks 
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are large, and they agree with quantum calculations at 
very high and very low frequencies, but more generally 
the heights of the Floquet peaks in the quantum calcu- 
lations remain mysterious. They will be explained using 
a semiclassical method described in the next section. 



IV. SEMICLASSICAL DESCRIPTION 

It is a general principle of quantum mechanics [25| that 
when in classical mechanics we add probabilities associ- 
ated with different paths leading to the same final state 
as in Eq. ()lf p . in quantum mechanics we add ampli- 
tudes. In the semiclassical approach, each amplitude is 
the square root of the classical density combined with a 
phase. In the present case, Eq. ([TT|) is replaced by 



P^ciPf)^\^sc {Pf,tf) 



with 



(13) 



X exp {i \Sj (p/, i/) jh - /ij7r/2 



-1/2 



(14) 



where we are again using = p{xo,tf). F{x()) is the 
envelope of the initial wave packet, either Fg{xo) in Eq. 
(ID) or Ftf{xo) in Eq. ([5]), and XQ{pf,tf) has the same 
meaning as in the paragraph above Eq. (jlip : trajectories 
that arrive at any one pf began from a large number of 
discrete XQ{pf,tf). 

Re-examining Fig.[3](a), and thinking about xo(p/,i/) 
as a smooth but multivalued function oipf, we divide the 
points XQ{pf,tf) into intracycle and intercycle groups, 
where a cycle is one period oipf- In Fig. [Slja) , we may say 
that the pair of points (a, /3) belongs to one cycle, the pair 
(7, 5) to another cycle, etc. Alternatively, we may say 
that the pair (/3, 7) belongs to one cycle, {5, e) to the next, 
etc. Summing over all the points a;Q {pf,tf) then means 
summing over points on distinct branches of xq (pf^tf) 
within a cycle, and then summing over cycles. Thus the 
index j may become a composite index, j — (b, c) where 
b is an integer labeling a branch within a cycle, and c is 
an integer labeling the cycle. 

J^{pf,tf) is a Jacobian, which in the present case is 
the same derivative defined in Eq. (|lip . 



^3 (Pf^tf) 



dpf {xo,tf) 



(15) 



Since pf \& a, periodic function of xq , the values of this 
derivative depend on the branches within a cycle, but 
do not depend on which cycle is examined: J^(b,c) iPf^tf) 
depends on the branch b but is independent of the cycle c. 
In Fig pa), {pf,tj) = J., {pf,tf) = ipf,tf) = 
while {Pf,tf) = Js {pf,tf) = Jq {pf,tf) = ... 

Sj {pf,tf) is a classical momentum-space action inte- 
grated along the path from XQ{pf,tf) to the final point. 




FIG. 7: (Color Online) Probability distributions of final mo- 
menta. The sharp peaks (blue online) are obtained by sum- 
ming over all branches of all cycles. Their heights are all 
multiplied by the same constant so that they are comparable 
to the other two curves. The oscillating curves are obtained 
by combining two branches of a single cycle, but with dif- 
ferent definitions of the cycle. The solid curve (red online) 
corresponds to a cycle spanning branches (6, c) in Fig. ETa), 
and the dashed curve (black online) is for a cycle spanning 
branches (c, d). Where those two curves intersect, the differ- 
ent cycles add in phase with each other, producing the sharp 
peaks. 



This integral is 

4" iPf^tf) 



xdp 



Edt 



X (xo, t) ^ ^ . ° ' ^ dt 



dt 



(16) 



E (t) dt 



There is a simple relationship between the values of 
^{b,c) (Pfi^f) for different cycles at fixed pf. Let label c 



increase with decreasing x^''^'; i.e., it increases by 1 with 
each successive cycle of the oscillating barrier. Then 



S{b,c+N) iPf,tf) = 5(b,c) {Pf, tf) + NAET, 



(17) 



where T is the period of one oscillation, N is the number 
of periods separating the cycles, and AE is the change 
of energy of the particle 



AE={p}^pI) /2m. 



(18) 



Finally, we introduce the Maslov index associated 
with each branch of a;o {pf^tf). The rule for determining 
it is given in Appendix B. Here let it suffice to say that 
in Fig. [3l^a), jlj can be taken to equal one on branches 
a, c, e, 5, ... and equal to zero on branches &, d, /, .... 

In our calculations, we compute the final momentum as 
a function of initial position pf {xQ,tf) = p{xQ,tf), then 
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for each pf we identify initial points Xq''^^ iPfytf) for 
branches b within a single cycle c. For each of them 
we find JbiPf,tf), Ab> and S(^i,,c) (Pf^tf) for that par- 
ticular cycle. We then calculate 5(6. c) {Pfjif) for other 
cycles using Eq. (1171) . and then compute the sum Eq. 
(fn)l numerically. Steps are also taken to correct the 
semiclassical approximation near divergent points, and 
the calculation is extended into the classically-forbidden 
regions; this procedure incorporates diffraction, or mo- 
mentum space tunneling, into the semiclassical dynam- 
ics. Derivation and additional details of the semiclassical 
method are given in Appendix [A] 

Terms in the sum over cycles add with incommensurate 
phases, and tend to cancel unless AE — 2ttK where K is 
any integer. This condition explains the Floquet picture 
introduced earlier: the momentum distribution becomes a 
"comb " function, with the "teeth " occuring at momenta 
that satisfy the commensurate phase condition, 



2m 



2m 



T 



(19) 



In Fig. [7| we show the absolute squares of two single- 
cycle wavefunctions, one using branches (&, c) (solid 
curve, red online) in Fig. [3Ka), the other using branches 
(c, d) (dashed curve, black online). These single-cycle 
probabilities intersect at momenta satisfying Eqs. © 
and (flQl) . The relative amplitudes of these intersections 
are determined by both the classical densities and the 
differences in momentum-space action, Eq. (|16p . among 
the paths contributing to the wavefunction at each p f . 

Fig.[8]shows quantities that determine the phase differ- 
ences and interference for three trajectories ending with 
the same final momentum. Figs. [SJa) and MJ^) show 
the position and momentum, respectively, versus time. 
Both plots show that particles see a decrease in velocity 
(and momentum) as they approach the potential bar- 
rier. Figs. [5Jc) and [HJd) illustrate the differences in 
the momentum-space action, Eq. (|16p . The differences 
in areas under the curves determine the phase differ- 
ences between pairs of trajectories. Interference asso- 
ciated with phase differences related to E{t) for different 
cycles (Fig. [HKd)) produces Floquet peaks. Phase dif- 
ferences between pairs of trajectories in the same cycle 
(Figs.[5{c) andlSI^d)) give the interference that determines 
relative heights of Floquet peaks. 

When we sum over cycles, the resulting probability is 
sharply peaked at the locations where the single-cycle 
probabilities intersect (Fig. [7|), and the heights of the 
peaks correspond to the relative magnitude of the single- 
cycle probability at these locations. Finally, we have an 
explanation for the relative heights of the Floquet peaks. 



V. CASE STUDIES 

In this section, we study three separate scattering cases 
for identical barrier parameters but different incident mo- 
menta: pure transmission, pure reflection, and mixed 
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FIG. 8: (Color Online) Quantities that determine the phase 
evolution and interference of three trajectories ending with 
the same final momentum. The solid (blue online), dashed 
(red onhne), and dotted curves (black online) correspond to 
the trajectories associated with {pf,xo) = P, 7, and 5 in 
Fig. ISja), respectively. One may think of the (/3,7) trajec- 
tories as being from a single cycle, with the 5 trajectory one 
cycle ahead of the (3 trajectory, (a) Position versus time. 
Each trajectory shows a decrease in velocity as the barrier is 
initially encountered near x = Q. (b) Momentum versus time, 
(c) x{t)dp(t) / dt term in the momentum-space action (Eq. I16|) 
versus time, (d) Energy term in the momentum-space action 
(Eq. p6p ) versus time. 



transmission and reflection. We also compare the full 
quantum results with the predictions of the semiclassi- 
cal approach and find relatively good agreement. While 
there is a large range of possible scattering behaviors that 
can be studied by adjusting the five input parameters of 
our model, these three cases capture most of the essential 
physics. 

Pure Transmission: The initial Gaussian wave packet 
is centered aX x — —1500, with (3 — 300, with initial mo- 
mentum po ~ 1-8, and with barrier parameters A — O.b 
and w = 0.1. This is the same case that was shown earlier 
in Figs. IHd), [31 [T] This initial momentum corresponds 
to an energy higher than the maximum amplitude of the 
barrier. It takes more than fifteen barrier oscillations 
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FIG. 9: (Color Online) (a) Final momentum vs initial po- 
sition for the Lj = 0.1, po = 1.8 case, (b) Comparison of 
classical (plotted upwards, black online), semiclassical (plot- 
ted upwards, blue online) and quantum-mechanical (plotted 
downwards, red online) momentum distributions. The hori- 
zontal lines in the upper portion of the graph correspond to 
the heights of the quantum-mechanical peaks. 



for the packet to pass over the barrier. There are two 
branches per cycle, as shown in Fig.lHl^a). The classically- 
allowed momentum values range from p/ ~ 1.2506 to 
2.1411. 

A comparison of Pg^ (pf) (plotted upwards, blue 
online), Pq (Pf) (plotted downwards, red online) and 
^ciPf) (plotted upwards, black online), is shown in 
Fig. ini^b). The semiclassical and quantum- mechanical 
results can be seen to agree well. The final probabil- 
ity has fifteen peaks within the classical envelope. Both 
the classical density and interference contribute to the 
relative heights of peaks. At least two non-negligible 
classically-forbidden peaks can be seen for momentum 
values on either side of the classical envelope. The semi- 
classical calculation has corrected divergent peaks near 
momentum turning points by using Airy forms of local 
wavefunctions (see Appendix E)) . 

Pure Reflection: We employ the same barrier parame- 
ters as in the previous case, but use an incident momen- 
tum of Po — 1.0, which corresponds to an energy equal 
to the minimum amplitude of the barrier. The barrier 
undergoes more than twenty-eight oscillations during the 
time the wave packet is interacting with it. There are two 
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FIG. 10: (Color Online) Same as Fig.[9]for po = 1.0. 



branches per cycle, shown in Fig.fTOla). with the classical 
envelope ranging from pf « —1.5043 to —0.6825. 

A comparison of P^q (pf) (plotted upwards, blue 
online), Pq (Pf) (plotted downwards, red online) and 
Pcipf) (plotted upwards, black online), is shown in 
Fig. [Toy b) , again with good agreement between the semi- 
classical and quantum-mechanical results. The final- 
momentum probability has nine peaks within the clas- 
sical envelope. We see at least three non-negligible peaks 
for classically-forbidden momentum values less than the 
minimum of the classical envelope, but only one non- 
negligible peak for forbidden momentum values greater 
than the maximum value of the classical envelope. This 
is because peaks are more closely spaced for large abso- 
lute momenta than for small absolute momenta, because 
they are equally spaced in energy. The exponential de- 
cay of the wavefunction again makes the peaks negligible 
outside the region shown. 

Mixed Reflection and Transmission: We implement 
the same barrier parameters as in the previous cases, but 
use an incident momentum of pq = 1.4142, which cor- 
responds to an energy between the minimum and max- 
imum of the barrier amplitude range. In this case, the 
wavepacket is partially reflected and partially transmit- 
ted. The periodic relationship between final momen- 
tum and initial position is more complicated in this case. 
Figs.fTlTa) and fT^ a) show the reflected and transmitted 
portions of the trajectory ensemble, respectively. Some 
classically-allowed final momenta have as many as six 
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FIG. 11: (Color Online) (a) Reflected portion of final mo- 
mentum vs. initial position for the u} = 0.1, po = 1.4142 case, 
(b) Semiclassical (plotted upwards, blue online), quantum- 
mechanical (plotted downwards, red online), and classical 
(plotted upwards, black online) final-momentum probabilities 
for the reflected portion of the wavepacket. 



FIG. 12: (Color Online) (a) Transmitted portion of flnal mo- 
mentum vs. initial position for the oj — 0.1, po = 1.4142 case, 
(b) Semiclassical (plotted upwards, blue online), quantum- 
mechanical (plotted downwards, red online), and classical 
(plotted upwards, black online) flnal-momentum probabilities 
for the transmitted portion of the wavepacket. 



interfering trajectories within each cycle. The classical 
envelope ranges from p/ ~ —1.6730 to 1.8987. 

Comparisons of the reflected and transmitted portions 
(pf) (plotted upwards, blue online), Pq {pf) (plot- 
ted downwards, red online) and P^{pf) (plotted up- 
wards, black online) are shown in Figs. fTTlb) and [T^ b). 
respectively. Every extremum in the pf{xQ,tf) graph 
gives a "turning point'' or caustic, at which Pc (pf) di- 
verges. The classical amplitude is markedly higher for 
larger momentum values in both the reflected and trans- 
mitted portions of the wave packet; consequently, the 
semiclassical and quantum-mechanical final-momentum 
distributions have their largest peaks in these regions. 
Agreement between semiclassical and quantum methods 
is less precise in this case, particularly where turning 
points are close together. 



VI. PROPOSED EXPERIMENT 

The theoretical predictions of the previous sections can 
be tested experimentally with the macroscopic wavefunc- 
tion of a BEC serving as the atomic wave packet. While 
the BEC does not have to be strictly ID, the use of a 



highly elongated BEC, confined in an optical dipole trap, 
simplifies the experiment. Furthermore, the BEC should 
be non-interacting since collisions between particles are 
not included in our calculations. A non-interacting BEC 
can be produced by employing a magnetic Feshbach res- 
onance. A number of alkali atoms, such as ^^Rb [2^ 
and ^^K [131, have been cooled to quantum degeneracy 
and also have a Feshbach "zero", a magnetic field which 
produces a null scattering length due to a nearby Fesh- 
bach resonance. In the non-interacting limit the Thomas- 
Fermi approximation no longer applies, and a harmoni- 
cally confined BEC has the Gaussian wavefunction of the 
trap ground state. 

The elongated BEC provides a wide, quasi-lD, Gaus- 
sian wavepacket, while a tightly focused blue-detuned 
laser serves as an optical dipole barrier with Gaussian 
shape, and with amplitude proportional to the laser in- 
tensity. Instead of launching the atoms towards the bar- 
rier, the opposite is more convenient in an experiment: 
The dipole barrier is swept through the stationary BEC, 
so that in the reference frame of the moving barrier the 
theoretical treatment still applies. Figure [T2] shows an 
optical circuit for generating a translating laser barrier 
using an acousto-optic modulator (AOM): the amplitude 
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FIG. 13: (Color Online) Proposed experimental implemen- 
tation of oscillating barrier scattering. An AOM generates 
a barrier laser beam whose deflection angle is controlled by 
the RF drive frequency. The beam rotation is converted to 
a translation (grey arrows) by a lens, which also focuses the 
beam to produce a narrow dipole barrier. A combination of 
lenses then inserts the laser barrier onto the BEG (atom sym- 
bols), ki and ^2, and associated arrows, indicate the wavevec- 
tors for the Bragg-Raman spectroscopy laser beams. The in- 
set shows a sample momentum and velocity distribution for 
the proposed Bragg-Raman spectroscopy experiment with a 
^^K BEG. 



and frequency of the radio-frequency (RF) drive for the 
AOM control the amplitude and position, respectively, of 
the barrier. 

The momentum components of the reflected and trans- 
mitted BEC wavefunction will generally be too close to- 
gether to be distinguishable by time-of-flight imaging. 
Instead, Bragg spectroscopy [28| can be used to measure 
the velocity distribution of the scattered atomic packet. 
Bragg spectroscopy is performed by briefly shining two 
laser beams on the scattered atoms, as shown in Fig. |T31 
When the two lasers are detuned from each other by a 
frequency S = Sq + {ki — ■ v, where fci and k2 are 
the wavevectors of the two incident Bragg beams, then 
only atoms with velocity v are given a two-photon mo- 



TABLE I: Table of proposed experiment parameters. 



Parameter 



Value 



atomic state 
Feshbach zero 
BEC width P 
BEC velocity 
barrier width a 
barrier amplitude Uo 
barrier mod. ampl. A 
barrier mod. frequency to 



\F = l,mF = +1) state of ^^K 
350 G 
10 fim 
12.9 mm/s 
2.5 nm 
197 nK 
1 

2n X lAkHz 



mentum kick h(ki — k^)- The energy imparted to an 
atom at rest by two-photon recoils determines the base 
detuning = hk'^ /{2m). The kicked atoms are spectro- 
scopically tagged with a Raman process [29] that changes 
their hyperfine level, accomplished by adding the hyper- 
fine ground level splitting to the base detuning Sq. The 
Raman-selected atoms are detected by absorption or flu- 
orescence imaging on the D2 line cycling transition. 

We summarize the main parameters of the proposed 
experimental implementation in table HI We consider a 
BEC of '^^K atoms in the \F ~ 1,tof = -1-1) hyperfine 
ground state, which has a vanishing s-wave scattering 
length at 350 G [23|- ^ red-detuned optical dipole trap 
produced by a 1 W 1064 nm laser focused to a 1/e^ diam- 
eter of about 120 /im will confine the BEC with a Gaus- 
sian density profile and an axial width of /3 = 40 = lO/im. 
A blue-detuned Gaussian barrier can be produced with 
a 532 nm laser focused to a radius of cr = 10 = 2.5/im 
(waist radius of 5 fim) with a barrier amplitude of Uq = 
1 = 197 nK. Translating this barrier at a velocity of 
12.9 mm/s (corresponding to an incident momentum of 
Po = 2 for particles of mass ra — 1 = 6.5 x 10^^^ kg), 
while modulating it at uj^arrier — 0.35 = 27rxl.4 kHz 
with a modulation strength of A=l, produces a purely 
transmitted wavepacket with the final momentum distri- 
bution shown in the inset of Fig. [131 The velocity peaks 
of the distribution have a half width at half maximum 
of Av wO.l mm/s, determined by the axial extent of the 
BEC. This velocity spread requires a laser frequency dif- 
ference stability on the order of 2Av/X «250 Hz (A = 767 
nm for ^^K), which is within the practical resolution of 
Bragg spectrocopy [s^]. Furthermore, we note that the 
axial confinement of the BEC does not play a significant 
role, since the trap has an axial oscillation frequency of 
f axial ~1 Hz, which is considerably slower than the time 
scale of the scattering process. 



VII. CONCLUSION 

In summary, we have studied scattering from an 
amplitude-modulated Gaussian barrier, and determined 
the final momentum-space probability distributions us- 
ing classical, semiclassical, and quantum formalisms. We 
find that classical mechanics defines the boundaries of 
a classically-allowed region of final momenta. Quan- 
tum calculations show: (i) the probability that particles 
end with momentum outside the classically-allowed re- 
gion is small; (ii) the momentum distribution is peaked 
at momenta consistent with Floquet's theorem; (iii) the 
heights of the Floquet peaks vary widely and seemingly 
erratically. Semiclassical calculations show that (a) for 
any final momentum inside the classically-allowed region, 
many classical paths arrive; (b) interference of waves 
propagating along these paths produces peaks consis- 
tent with Floquet theory, and determines their heights. 
Specifically, inter-cycle interference leads to discrete final 
momentum states, while intra-cycle interference deter- 
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mines the peak heights. Finally, momentum-space tun- 
neling leads to difFractive population of momenta beyond 
the classically-allowed bounds. 

The semiclassical and full quantum propagation for- 
malisms employed in this work are well suited for study- 
ing scattering from a turnstile pumping potential formed 
from two separated barriers, amplitude-modulated 7r/2 
out of phase from each other. Classically, such a poten- 
tial displays strong signatures of chaos, with quantum dy- 
namics well suited to the type of semiclassical treatment 
developed in this paper. Such a treatment is essential 
for understanding the quantum and classical aspects of 
particle pumping in a turnstile pump, since interference 
and tunneling can be selectively included. Moreover, the 
scattering theories developed in this work can also be 
extended to examine spatial tunneling through narrower 
barriers, and scattering from a potential well. 
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Appendix A: semiclassical Analysis 

We give here details and derivation of the semiclassical 
formulas used in Sec. IIVI Most of the theory is similar to 
methods we have used in earlier papers 3lJ-|41| , but some 
aspects of the present system are different. In most of our 
earlier work, we have studied stationary fixed-energy sys- 
tems; only [35l438| dealt with time-dependent potentials. 
In the present case, the initial and final conditions are, 
from semiclassical perspectives, a little unusual. At the 
final time, we want a semiclassical approximation in mo- 
mentum space. However, at the initial time, we cannot 
use a semiclassical approximation in momentum space, 
though we can in configuration space. Furthermore, the 
sum over cycles of the oscillating barrier is different from 
previous work. 



1. Local Wavefunction 

Recall that we have an oscillating Gaussian barrier 
with a wave packet approaching from the left. At an 
initial time to, the wavefunction for x <C (far to the left 
of the barrier), is given by 



*o(a;,io) = F(a;)e'(P°^-^°*°)/'', 



(Al) 



where F{x) is a function describing the envelope of the 
initial packet in (x, t) space. We include time and en- 
ergy as canonical variables, expanding the phase space 



for the system. For reasons that will become clear, one 
regards t as a canonical momentum, and E as a canonical 
coordinate, q — (x, E) and p — {p, t). 

Then defining an effective Hamiltonian, , given by 



P 

2m 



+ U{x,t)-E, 



the equations of motion are 



(A2) 

(A3a) 
(A3b) 
(A3c) 
(A3d) 
(A3e) 
(A3f) 



where t is a "timelike" progress variable along the trajec- 
tories, and is related to t in the Schrodinger Equation via 
T = to + 1. We call S the classical action along the tra- 
jectory, and S can be thought of as a "momentum-space 
action" along the trajectory. The form of equations (jA3cP 
and (jA3dl) justify the indentification of as a canonical 
coordinate and t as its conjugate momentum. 

We want to compute the probability that the particles 
end with a given final momentum, using the momentum- 
space wavefunction 5" (p, t). Therefore, we want a semi- 
classical approximation in momentum space. However, 
since we have chosen an initial distribution with very 
small momentum spread, the initial wavefunction in mo- 
mentum space is nearly a delta function, which cannot 
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FIG. 15: (a) Slice of Lagrangian manifold at small time, (b) 
Periodic final momentum as a function of initial position, (c) 
Final momentum, pf — p (xq, Tf), as a function of final posi- 
tion, Xf = X {xo, Tf). This corresponds to the final-time slice 
of the Lagrangian manifold. 
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FIG. 16: (Color Online) Slice of Lagrangian manifold at an 
intermediate time. The numbers correspond to intermediate- 
time slices of different momentum charts, which are separated 
by local extrema in the function p — .^{Xji) for fixed time, 
denoted by large circles. For every given momentum (e.g., the 
dashed line), there are many corresponding values of x. 



be described by a semiclassical approximation. There- 
fore, in order to calculate the desired momentum-space 
wavefunction, we start our calculation in (z, t) space, and 
later transform to (p, t) space. 

The first step in constructing a semiclassical wavefunc- 
tion is is to propagate trajectories from a line of initial 
conditions. We choose the line of initial conditions to 
have a constant starting time = 0, variable starting 
position X covering the domain of the initial packet, and 
a fixed initial momentum pg- The resulting trajectories 
sweep out a two-dimensional surface called a Lagrangian 
manifold in the four-dimensional (x,p,E,t) phase space. 
A typical Lagrangian manifold for this system is shown 
in Fig. [H 

Integration of trajectories with respect to t gives a 
relationship between {xq,t) and {z,t), where z is any 
dynamical variable x,p,E, S, or S. From our choice of 
to — 0, t is simply equal to t, and x is the point at which 
the trajectory arrives at time t ~ t. We may think of 
each of these quantities as a function of the initial variable 
Xq and the progress variable r, e.g., x{xo,t),p{xo, r), etc. 

We define a Jacobian, 



J{xo, t) = det 



dix,t) 
d{xo,T) 



dx 
dxo 



7^ (A4) 



with Jo = J {xo,0) = 1. This Jacobian is a single- valued 
function of {xq , r) . For r not too large (and x not too 
far from xq) there is an invertible relationship between 
(xqjt) and (x,i); i.e., we may consider (xo,t) as a func- 
tion of {x,t). With this relationship, we may also con- 
sider the position-space action S and Jacobian J to be 
functions of {x,t), 



S{xo,t) ^S{xoix,t),T{x,t)) ^S{x,t) 
J{xq,t) ^J{xaix, t), t{x, t)) =Jix, t) (A5) 

We may use these functions in the primitive semiclassical 
approximation for the [x, t) space wavefunction 



^sc(x,t) = ^'o(xo,T = 0) 



Jo 



J{x,t) 



1/2 



(A6) 



where (xq, t) are considered to be functions of (x, t). The 
initial Maslov index has been set equal to zero, and 



vE'o(xo,r = 0) = F(xo)e'f°-°/^ 



(A7) 



where (xo,t) are again considered as functions of {x,t). 

As the trajectories are propagated forward in t, they 
come to the barrier region, where p is no longer con- 
stant, and we may use (p, t) locally as independent vari- 
ables to describe the Lagrangian manifold, as shown in 
Figs.[ll[Tii;a) andlTSjc), and Fig. [TH 

A "momentum chart" is a region of the Lagrangian 
manifold that has a diffeomorphic projection to momen- 
tum space, (p, i). In Fig. [Tni a constant-time slice of 
the Lagrangian manifold is shown. For each value of p, 
there are many corresponding values of x; each can be re- 
garded as a "branch" of a multivalued function, and each 
is a constant-time slice of a momentum chart. 

We transform to the momentum-space wavefunction 



(A8) 
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We evaluate the integral for the part of the wavefunction 
that corresponds to the initial momentum chart by using 
the stationary phase approximation. We use the func- 
tion .'^{x,t) = dS/dx to describe the Lagrangian man- 
ifold, and p is the independent variable in vE'(p, t). When 
we substitute the semiclassical approximation (jA6|) into 
(|A8|) . each classically-allowed p has a stationary phase 
point, X, where p — .^{x,t), i.e., where the line p — 
constant intersects the Lagrangian manifold, as shown in 
Fig. [in] for p = 1.7. In evaluating the integral, we also 
make use of the momentum-space action, defined in Eq. 
(IA3fp . and define a momentum-space Jacobian 



on this slice of the manifold, they can be connected by 
a path on this slice. At each point that the path passes 
through a momentum turning point, the Maslov index 
changes by ± 1 , and we use the following rule to determine 
the increment. This rule applies if the {x,p) plane is 
drawn in the most usual way, with x increasing to the 
right and p increasing upward. When the path passes 
through a momentum turning point that separates the 
momentum chart from the j*'' momentum chart, then 

Mi = Mi + I7 if path curves right (CW) (Al3a) 
jlj = if the path curves left (CCW) (A13b) 



■th 



J{xo,t) = det 



dip,t) \ ^ dp_ 
d{xo,T)J dxo 



(A9) 



The locally invertible relationship between (p, t) and 
(xo,r) allows us to consider S{xq,t) and J{xq,t) to be 
functions of {p,t), i.e., 



S{xo{p,t),T{p,t))^S{p,t) 

J{xo{p,t),T{p,t)) ^ J{p,t). 



(AlO) 



With these definitions, the stationary phase approxima- 
tion in the initial momentum chart yields 



-1/2 



(All) 



Generally, for every momentum chart of the La- 
grangian manifold, there is a comparable term contribut- 
ing to the momentum-space wavefunction. We write the 
local, primitive form of the momentum-space wavefunc- 
tion for each momentum chart as 



^,{p,t) = F{xl{p,t)) 



1 



1/2 



X exp 



iSj{p,t) injlj 



(A12) 



where flj is the Maslov index for the given momentum 
chart. 



a. Maslov Index 

Here, we state the rule for the Maslov index for each 
momentum chart. As indicated in Fig. 1161 momentum 
charts are separated by momentum turning points, which 
are extrema of locally-defined functions p — -^{x, t) for 
fixed i, i.e., points where d.'^{x,t)/dx — 0. 

Each time any path on the Lagrangian manifold passes 
through a momentum turning point, the Maslov index 
changes. In Fig. ITsy c). we show a slice of the Lagrangian 
manifold at the final time tf. If we take any two points 



where CW and CCW denote clockwise and counter- 
clockwise, respectively. 

The exp (— i7r/2) term in the primitive wavefunction for 
the first momentum chart, (|All| . corresponds to jlj = 1 
in (jA12p . All other Maslov indices for the remaining 
momentum charts are constructed relative to it, using 
(|AT3a| and (|AT3bl) . 

For the two paths shown in Figs. fTSl c) and [TBI moving 
from left to right, the Maslov index increases at every 
maximum, and decreases at every minimum. 



h. Corrections Near Momentum Turning Points 

The primitive semiclassical wavefunction diverges at 
momentum turning points, where where J7b(p, t) vanishes. 
To correct this, we construct an alternative way of writ- 
ing the primitive wavefunction, which will be valid near 
momentum turning points in classically-allowed regions. 
We then match this form of the wavefunction to the Airy 
function and its derivative, in order to extend the semi- 
classical approximation into classically-forbidden regions 

sa. 

We start by adding the primitive forms of the wave- 
function, (jA12p . for two successive momentum charts, 
and we denote this wavefunction \l/m+n(p, t)- We intro- 
duce the following notation 

A{p,t) = \j{p,t)'^^'^ (A14a) 

A>S(p, t)=Sn{p,t)^ S,n {p, t) (A14b) 

t) = [5„ t) + Sm {p, t)] /2 ( Al4c) 

AA{p, t) = A„{p,t)- A,n [p, t) (A14d) 

A(p, t) = [A„ (p, t) + A„ (p, t)]/2 ( Al4e) 

AF{xo{p,t)) = Fn{xo{p,t)) - Fmixoip,t)) (Al4f) 

¥{xo{p,t)) = [Fr,ixoip,t)) + Fm{xo{p,m /2, (A14g) 

where the m and n subscripts denote the momentum 
chart with the lower and higher Maslov index, respec- 
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tively. We use these definitions to write 
^,n+n(P,t) = 2exp — — 



X <^ AF + 



AAAF\ 



+ 



4 J 
fAAF AA¥ 



e-"/4 sin 



V 2 2 



6-*^'^/* COS 



ASip,t) 


TT 


2h 


^4 


AS{p,t) 


TT 


2h 


^4 



(A15) 



We match the separate terms of (jA15[) to the first-order 
asymptotic forms of tire Airy function and its derivative, 
respectively, so that we may write (jAlSp as 



*m+n(p, t)^C (p, t) Ai(-Z(p, t)) 

+ D{p,t) Ai'(-z(p,t)), 



(A16) 



where 



AAAFl 



D 



^-1/2 (z(p,t))-^/^ 

-2exp(^(|^%i:-f))[M£ + A^] 

2/3 



Ah 



(A17b) 
(Al7c) 



We use wavefunctions of the form of (|A16[) in 
the classically-allowed regions near momentum turning 
points, where (jA12[) is not valid. 



c. Classically-Forbidden Regions 

One can show that if the momentum turning points 
are quadratic maxima or minima, the following functions 
vary linearly with p near the turning point p 



[ASip,t)r^^ 


cx {p 


-p) 


(Al8a) 


E{p,t) + S{p,t) 


(X {p 


-p) 


(A18b) 


[A{p,t)]-' 


cx {p 


~p) 


(A18c) 


[AAip,t)f 


cx {p 


^p)- 


(AlSd) 



We continue these quantities into the classically- 
forbidden regions using these linear approximations. To 
obtain values for F{xq{p, t)) in these regions, we extrapo- 
late xq into the classically-forbidden regions, and use it to 
evaluate F{xo{p,t)). This extrapolation yields complex 
values of xq . 



2. Global Wavefunction 

We denote as "branches" the regions separated by mo- 
mentum turning points in p{xo,Tf)^ i.e., regions sepa- 
rated by points where dp{xo,Tf)/dxo = 0. We define 
a "cycle" as one barrier oscillation, i.e.. one period of 

PixQ,Tf). 

We want to construct a final wavefunction that is valid 
in both classically-allowed and classically-forbidden re- 
gions. We have seen that each momentum chart con- 
tributes a term to the final wavefunction, so our first 
step is to construct all local wavefunctions. 

We will illustrate the steps necessary to construct the 
final wavefunction for the simplest case, like that shown 
in Fig. [3]Ja), which contains two branches per cycle. We 
must determine the regions of validity of the two forms of 
the wavefunction, Eqs. (jA12p and (jA16l) . for all branches. 
Due to the periodicity of final momentum and initial po- 
sition, we can do this for a single cycle only, as Eqs. 
(|AT2| and ((A16| are valid in the same regions for the i* 
branch within every cycle. Further consequences of this 
periodicity are discussed in App. [B] 

We choose the cycle spanning branches (a, b, c) in 
Fig. EI a) . We start with branches a and b, and construct 
the primitive form of the wavefunction by adding Eq. 
(|A12p for the two branches. We then construct ^a+b{p, t) 
via Eq. (IA16p . These two forms of the wavefunction are 
valid in different but overlapping regions, and we com- 
pare the two to determine the region of validity for each. 
This comparison shows that the Airy form is valid in re- 
gions D and E in Fig.EJa) {p < 1.66). In region D (1.36 
< p < 1.66), where both forms of the wavefunction are 
valid, we use a switching function that varies between 
and 1 to weight each form, and use a linear combination 
of the two. We then use 



^a+b{p,t) =/i (p) [Airy form] 

+ (1 — /i (p)) [Prim, form] 



(A19) 



as the local wavefunction for branches a and b in regions 
C, D, and E, where /i is the switching function; /i — 
at the boundary between regions C and D. "Airy form" 
and "Prim, form" in Eq. (|A19|) refer to ^a+bip,t) calcu- 
lated via Eqs. (|A16p and (|A12p . respectively. 

We repeat this process for branches b and c, and find 
that the Airy form of this wavefunction, ^b+c(p, i)j 
valid in regions A and B {p > 1.77). Both forms of the 
wavefunction are valid in region B {1.77 < p < 2. 01). We 
use a switching function in region B to weight each form 
of the wavefunction, and use a linear combination of the 
two. We use primitive semiclassical wavefunctions for all 
branches in region C (1.66 < p < 1.77). 

With knowledge of where each branch's primitive and 
Airy forms of the local wavefunction may be used, one 
may construct a final wavefunction, which is a linear com- 
bination of all local wavefunctions. For cases with more 
than two branches per cycle, a more elaborate version of 
the same process is used. 
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Appendix B: Semiclassical Implications of 
Periodicity 

An initial wavefunction that is long in position space 
needs many oscillation cycles to pass through the barrier 
region. Semiclassically, this means that the summation 
of the primitive wavefunction ^j{p,t) (Ea. (|A12p ) over 
the momentum charts j involves a sum over trajectories 
with initial xq values extending over numerous oscillation 
cycles of p{xq, tj), as seen in Fig. [T5lb). This creates in- 
terference of trajectories belonging to different cycles of 
oscillation. This inter-cycle interference constructively 
enhances final momentum values satisfying AE = huj, 
consistent with Floquet theory. Here we clarify how this 
constraint arises semiclassically and derive explicit for- 
mulas for the resulting momentum-space wavefunctions. 
The following discussion refers to the classically-allowed 
regions, but its validity could be extended to include re- 
gions near turning points, and classically-forbidden re- 
gions, by using the appropriate Airy forms of local wave- 
functions. 

Let L denote the initial interval of a;- values over which 
the initial wave packet is defined. We consider all those 
trajectories ending with a given value oi pf and begin- 
ning with any initial xq in L. We further restrict atten- 
tion to trajectories whose final point x{xo,Tf) is suffi- 
ciently far outside the barrier region that the potential 
is essentially flat. This is appropriate when most of the 
wave packet has either reflected from or passed through 
the barrier region. Now, we choose some interval / of 
length poT/m, corresponding to one oscillation period 
T = 2tt/uj, within L. Label all those trajectories ending 
at pf which have xq inside / with an index b as above; 
i.e. the initial position for each such trajectory is labeled 
Each Xq is one member of an entire family of initial 



positions x, 
note that x, 



{b,c) 


(6,0) 



cpoT/m, indexed by an integer c; 
Xq. Thus, b (branch) labels trajectories 



within one oscillation cycle of Fig. [TST bl. and c (cycle) 
distinguishes trajectories between different oscillation cy- 
cles. 

The primitive form of the momentum-space wavefunc- 
tion is given by summing (|A12[) over the double index 
J = {b,c): 




1 1 


Gaussian 


1 N = 1 1 


^ rectangular ^ 



FIG. 17: Plot of |Z)|^ for an initial rectangular [lower curve, 
Eq. (|B7|l ] and Gaussian [upper curve, Eq. (IB10|) [ initial packet 
profiles. For the rectangular case, N = 3, 10, showing conver- 
gence to delta functions at integer values of e. The widths 
/3 of the Gaussian packets are chosen to match the standard 
deviations of the corresponding rectangular packets. 



index, i.e. 



J(b,c){P,t) = J(b,0)iP,t) = JbiP,t), 
P-{b,c) = A(b,0) = Mb- 



(B2) 



The actions too can be related to one another. Consid- 
ering first 5(h^o) a.nd S(^b,i) , the (6, 0) and (6, 1) trajectories 
follow the same path in the barrier region, but the (6, 1) 
trajectory spends one more cycle to the left of the bar- 
rier, whereas the (6, 0) trajectory spends one more cycle 
to the right. Hence by Eq. (jA3f|) 



S(b,i){p,t) - S(i,^Q){p,t) = AET, 



(B3) 



where AE = p§/2 is the energy gained (or lost) by 

the trajectory due to scattering from the barrier. Since 
AE does not depend on the indices b or c, we conclude 
that 



= 5, 



(6,0) 



cAET = Sb + cAET. 



(B4) 



Eqs. (|B2p and (IB4[) provide an efficient method for con- 
structing terms when computing the semiclassical wave- 
function. Rather than directly integrating the entire line 
of initial conditions L, one only needs to integrate tra- 
jectories for initial conditions within one cycle, e.g. the 



interval /, and construct S, 



{b.c) 



for other branches via 



B4l). The semiclassical sum can thus be rewritten as 



b c— — oo 



X exp 



J{b,c){P,t) 

iS(b,c){p,t) i'^fj'{b,c)\ 



1/2 



(Bl) 



^{p,t) = 

J2^b{p,t) 



b 

where 



Jbip,t) 



/ iSb{p,t) infib] 
exp ; — ) (B5) 



Here, we allow c to range over all integers, since the initial 
profile F{xq) serves to effectively eliminate any trajecto- 
ries that begin outside L. Since two trajectories with 



initial positions Xq''"'^'' 



and X, 



(b,c') 




having the same b in- 



dex, differ only in their (uniform) motion outside of the 
barrier region, they have the same Jacobian and Maslov 



Db{p,t)= J2 i^(4(P,^)-cpoT/m)e'=^^^/^ (B6) 



In most of our calculations, we perform this sum numer- 
ically. However, in some cases, the sum can be expressed 
in closed form. 
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FIG. 18: (Color Online) For the "elevator" model, the change 
in energy vs. phase of the oscillation, (j) (solid line, blue on- 
line); potential energy V{x = Q,t = Q) (dotted line, black 
online); p (dash-dot line, red online), and tb (dashed line, 
green online). Parameters are po = 2.0, Uo — 1, A = 0.5, 
m = l,a; = l, L = 0.1. The energy change is plotted as 
AE/{mu}L/po). The points on the change in energy curve, 
going from left to right, correspond to Eq. (|C2c|) evaluated 
with our chosen parameters for (p = {(f)y,<j)<) (see Eqs. (|C3|) 
and (|C4|| '). respectively. 



We consider Eq. (|B6[) for two initial packet profiles, 
rectangular and Gaussian. Considering the rectangular 
profile first, take F{xo) = Fq constant over an interval of 
length NpoT, corresponding to N oscillation cycles, and 
F{xo) = outside this interval. Then Eq. (jB6[) can be 
rewritten as 



D,ip,t) = Dip) 



Fo e^^^*" 



c=0 



sin(7reiV) 
sin(7re) 



(B7) 



where e = AE/{hijj). Since D{p) does not depend on b, 
Eq. (jBSP factors into the product oi D{p), involving only 
a c sum, and a quantity involving only a sum over b. 



^{p,t) 
D{p)^ 



b 



Mp,t) 



1/2 



exp 



iSb{p,t) in fib 



(B8) 



As the length of the initial wave packet goes to infinity 
(i.e. N goes to infinity), approaches a comb of delta 
functions according to 



lim y Sie-k). 



(B9) 



Thus, the scattered wavefunction obeys AE = ktku, in 
agreement with Floquet theory. Convergence to the delta 



FIG. 19: (Color Online) Energy changes vs. phase for the 
same parameters as in Fig. 1181 except L — 0.3, 1.0, 3.0 and 
10. 



functions is illustrated by the lower curves in Fig. [171 
which show D{p) (Eq. (jB7|) ) as a function of e for = 3 
and 10. 

Considering the Gaussian profile next, we now take 
F{xo) equal to Fq (xq) in Eq. Then in the limit of 

a long packet {/3 >> poT), Eq. (|B6|) reduces to 

Dbip,t) = Dip) = -,=±^e,{.e,, 
where 6*3(2:, q) is a Jacobi theta function [43 



(BIO) 



93{z,q)^l + 2j2 9" cos(2nz) 



(Bll) 



The upper curves in Fig. [TTI illustrate D (Eq. (jBlOll ) as 
a function of e. As in the case of a rectangular initial 
condition, D converges to a comb of delta functions as 
the initial packet width increases. Unlike the rectangular 
case, however, there are no higher order peaks between 
the primary peaks at integer values of e. This agrees 
with the results presented in the paper (Figs. [ WT^ . 
which also show no higher-order peaks between the pri- 
mary Floquet peaks. 



Appendix C: Boundaries of Classically-Allowed 
Regions 

It would be nice to obtain some simple estimates of 
the maximum and minimum classically-allowed energy 
change. This turns out not to be as easy as we might 
wish. The simplest model is an "elevator:" 



Viix,t) 



Uo{l + Asm{Lut + (f>)) ,0<x<L 







otherwise 



(CI) 
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A particle of mass m and initial momentum po arrives at 
X — at time t = 0. If at that instant its kinetic energy 
is less than 1^1(0,0), then the particle is reflected with 
momentum — po- Otherwise it hops onto the elevator, 

1/2 

traverses it with momentum p' = [pq — 277114(0,0)] , 
and arrives at the end of the elevator at time tb = mL/p'. 
There it hops off, gaining potential energy Vi(L,ifc), so 
the final energy and the change in energy are 



Ef = 



2m 



Vi{L,tb) 



AE = Vi{L,tt)~ViiO,0) 

= AUq [sin (wtf, + — sin ( 



(C2a) 

(C2b) 
(C2c) 



The maximum possible range of AE is ±2 AUq. It is also 
important to note that tb depends on 0. 

Intuitively we expect that if ujtb is small, then the par- 
ticle will gain the most energy if it arrives at the barrier 
when the elevator is most rapidly rising, i.e., if cf) — 0. 
This is a respectable guess; however, if it arrives a bit 
later, then it will spend a longer time on the elevator, 
and thereby gain more energy. Likewise, we may expect 
that it will lose the most energy if it arrives when the ele- 
vator is falling most rapidly, (p — n. However, if it arrives 
a bit earlier, then again it stays longer on the elevator, 
and so it loses more energy. 

A graph of AE vs. (j) is shown in Figs. [18] and [19] for 



small Lotb- The maximum increase in energy occurs when 
= 0> , where 



mAUo 



(C3) 



and the greatest decrease occurs when <f> = <f>^. where 



(/)< « TT 



{pI - 2mUn) 



(C4) 



The change in energy predicted by these values of are 
shown in Fig. [TS] For wider barriers, the behavior be- 
comes more complex. 

There are also other solvable models, such as 



V2 {x, t) 



Vo{t)~\x\ ,Q<\x\<Vo{t) 







otherwise 



(C5) 



and 



V:i{x,t) 



where Vq (0 
complicated. 



(C6) 



, otherwise 

C/q (1 + ^ sin (cjt)), but they are more 
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